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1 Introduction 

Integrable models with quantum superalgebra symmetries have been the focus of recent 
studies JI], 0, 0, |, H, in the context of strongly correlated fermion systems, a subject 
of high-profile international research activity because of their relevance to high- Tc super- 
conductivity. The investigations to these models have largely been carried out within the 
framework of QISM and Bethe ansatz method. The exception are the works in [0, where 
the algebraic analysis method, developed in [R P| and generalized in O, O, [TB|, IH, O, 



was used to diagonalize the supersymmetric t-J model and its multi-component version 
directly on an infinite lattice. 

The algebraic analysis method & H], which we will call the vertex operator method. 



was formulated with the help of the level-one q-vertex operators |]10| and highest weight 



representations of quantum affine algebras. The vertex operator method was later ex- 



tended in |T6| to treat integrable models with boundary interactions ||I7| , |lq] . It was 
shown in how the space of states of the boundary XXZ spin-| chain on a semi-infinite 
lattice can be described in terms of level-one q-vertex operators of Ug{sl2), and how the 
correlation functions can be computed by the vertex operators. Several other models have 



been analysed by means of this approach |2^, ^ . 

In this paper, we study the q-deformed supersymmetric t-J model with an integrable 
boundary. We will work directly on a semi-infinite lattice. As is known, the q-deformed 
supersymmetric t-J model on an infinite lattice (i.e. without a boundary) has as its 
symmetry algebra the quantum affine superalgebra [/g[sZ(2|l)] [0. On a finite lattice 
with diagonal boundary refiection K-matrices this model was solved in by the Bethe 
ansatz method. Here we adapt the vertex operator method. We will diagonalize the 
boundary model Hamiltonian directly on the semi-infinite lattice, and moreover compute 
the correlation functions of the boundary model. 

This paper is organized as follows. In section 2, we describe the vertex operator 
approach to the q-deformed supersymmetric t-J model on the semi-infinite lattice. In 
section 3, we study the bosonic realization of the boundary states associated with the level- 
one highest weight representation of [/g[s/(2|l)]. In section 4, we compute the correlation 
functions of the local operators (including the spin operator Sf ) and derive the difference 
equations which they satisfy. In appendix A, we review the bosonization of f/g[s/(2|l)] at 
level-one and the associated vertex operators. 



2 Boundary q-deformed supersymmetric t-J model 

2.1 q-deformed supersymmetric t-J model on a finite lattice 

In this section, we recall some facts about the q-deformed supersymmetric t-J model on 
a finite lattice. Throughout this paper, we fix q such that |g| < 1. 

Let V be the 3-dimensional graded vector space and Eij be the 3x3 matrix whose 
(i, j)-element is unity and zero otherwise. The grading of the basis vectors vi, V2, fa of 
V is chosen to be [vi] = [V2] = 1, [fa] = 0. Let V* be the dual space and {fiifo)"^-!} 
the dual basis vectors. Denote by Vz (resp. V^*'^) the 3-dimensional level-0 representation 
(resp. dual representation) of ?7g[sZ(2|l)] associated with V. Let R{z) G End{V ® V) be 
the R-matrix of t/g[s/(2|l)] with matrix elements defined by 

R{z){Vi (g) Vj) = J2 Ru{z)Vk ® Vi, Vf j, Vj, Vk, Vi G V, 
k,l 
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otherwise. 

The R-matrix satisfies the graded Yang-Baxter equation(YBE) onV ®V ®V 

Rl2iz)Rl3izw)R23iw) = R23iw)Ri3{zw)Ri2iz), 

and moreover enjoys: (i) initial condition, -R(l) = P with P being the graded permutation 
operator; (ii) unitarity condition, -Ri2(^)-R2i(7) = 1, where i?2i(^) = PRi2iz)P] and (iii) 
crossing-unitarity, 

sti 



R-^''''{z) (^{M ® l)R{zq-^){M ® l)y ' = 1 ® 1 



where 



2pi 



2p2 



2P3 



(2.1) 



The various supertranspositions of the R-matrix are given by 

\kl 
lij 



iT'^'izW = = R(z)%. 



Following Sklyanin , we construct the transfer matrix of an integrable finite chain, 
with an open boundary condition described by a refiection K-matrix K{z). Here K{z) is 
a solution of the graded refiection equation 



K2iz2)R2liZlZ2)Ki{zi)Ri2iZi/z2) = i?21 (^1 7^2)^1 (^1 )^12 (2:122)^2 (2:2) • 

With appropriate normalization, we can show that this K{z) obeys the relations 



(2.2) 



K(l) = l, 

K{z)K{z-^] 

Kiz)Kiz-'] 



(Boundary initial condition), 
1, (Boundary unitarity), 
1, (Boundary crossing — unitarity). 



(2.3) 



where K{z) is defined by 
K(z) = 



^i?(^2)g^_^)[i] + [i] + b-][/3] + [a][/3]^^(^-l^-l)^2p.^ 



(2.4) 
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The third relation is the graded extension of the boundary crossing -unitarity p|, |6H2|. 



The transfer matrix of the q-deformed supersymmetric t-J model on a finite chain 



with the open boundary condition is constructed from R{z) and K{z) via [|T^, Q 

Tl-{z) = stry,{K\z)r{z~')K{z)r{z)), (2.5) 

where K+{z) = K{-z-^q-^Y^M and 

r{z) = Rqi{z) ■ ■ ■ Rqn{z) e End{Vo ®Vi ■ ■ ■ (^Vn) 

is the double-row monodromy matrix. The supertrace is defined as str(A) = l)'*'^M• 
It can be verified that T^^{z) form a commuting family, [T^^{z),T^^{w)] = 0. The 
Hamiltonian of the boundary q-deformed supersymmetric t-J model is given by [|r^, |^ 

- Tz^^ ^'^^'=' - f^, + ~2d'z^^'^^'=' + K^) ' ^^-^^ 

where hj^j+i = Pjj+i£Rjj+i{z)\,=i. 

The transfer matrix ( p.5|) with diagonal reflection K-matrices was diagonalized by the 
Bethe ansatz method in [^. 

2.2 q-deformed supersymmetric t-J model on a semi- infinite lat- 
tice 

In this paper, we restrict ourselves to the diagonal reflection K-matrix of the form 



/ i^z 



K{z) = f{z) 



z~r 



V 




1 I , f{z) = 0(z|r) = (2.7) 

(p{z ^|r) 1 — rz 



where r is an arbitrary parameter which is related with the boundary interaction [0, {Wj . 
One can check that such a K-matrix satisfies the boundary unitarity and crossing-unitarity 

(PD- 

We now consider Hamiltonian (^.61) in the semi-infinite limit: 



i=i 

where A = |^/(-2)U=i acts formally on the left-infinite tensor product space 

■■■^V(g)V. (2.9) 

As mentioned in the introduction, the q-deformed supersymmetric t-J model on an 
infinite lattice has ?7g[s/(2|l)] as its symmetry algebra. Let V{fia) be the level-one irre- 
ducible highest weight f/g[s/(2|l)]-modules with highest weight /Iq, a G Z (see (|A.4|) and 



0). Consider the level-one vertex operators which are intertwining operators between 
V{fia) and V{fii3). It has been shown in that the following type I vertex operators $(z) 
exist, ^*{z) which interwine the level-one irreducible highest weight ?7q[s/(2|l)]-modules 



^z) : Vifi^) F(/ia-i) ® V;, $*(^) : \/(/i„) Vifia+i) ® Vf. (2.10) 
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(See Appendix A for more details about V{fia) and its associated vertex operators.) 
Therefore following |^ |TB|, we can write the transfer matrix of the q-deformed 
super symmetric t-J model on the semi- infinite lattice as 

= E q-''^^A^)Kl{z-\-'m{z-'), (2.11) 

where ^i{z) and ^j{z) are the components of the ?7g[s/(2|l)] vertex operators of type 
I (see ( |A.6|) ). We have used the exchange relations of vertex operators ( |A.14| ) and the 
definition of K{z) ( p.4| ) in the above equation. 

We remark that the transfer matrix Tb{z) given by (p.ll|) is an operator with the 
property 

T{z): y(/i„) ^ r(/ij, aeZ. 

The commutativity of the transfer matrix (|2.11|) , [Te(-2), Tb(w)] = , then follows from 
( |70^ ) and (^). Moreover by (|70^ ), (|70B| ) and ( [70BD , one can show 



Tb{1) = id, Tb{z)Tb{z-^) = id, (2.12) 
Tb{z)Tb{z-\-^) = id. (2.13) 

These relations correspond to the boundary initial condition, boundary unitarity and 



boundary crossing-unitarity (|2.3|) of the K-matrix, respectively. In terms of the transfer 
matrix, the q-deformed supersymmetric t-J model Hamiltonian on the semi-infinite lattice 
is given by 

H = ^Tb{z)1=i. (2.14) 
Following 1^, we define the local operators acting on the n-th site: 

= -«f*(l)«f,(l)(-l)[^], (2.15) 
4? = E(-l)^''^^^'^^^'"'^'"^'^m(l)4r'^'^-(l)' ^ = 2, 3, ■ ■ ■ • (2.16) 

m 

In particular, we have the spin operator SI 

SI = \{E^i^-E^j) = \ {$t(l)$,(l) -$;(1)$2(1)}. 

3 The boundary states 

In this section we construct the bosonic boundary state |a;r >b and its dual state b < 
r; a|, which satisfy 

TB{z)\a]r >B=\a]r >B, b < r]a\TB{z) = b < r]a\. (3.17) 

By ( |A.15|) and (|2.11| ), the above eigenvalue problem is equivalent to 

^i{z'^)\a-r>B = Y.Kl{z)^A^)\ci]r>B, (3.18) 

i 

B<r-a\^*{z){-l)^^ = Y.B<r;a\^*{z-^)Ki{z){-li'^. (3.19) 
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3.1 The boundary state in V{Ao) 

Firstly, we consider the boundary state |0 >bG V{fio) (or V{Aq)). As is shown in Ap- 
pendix A, V{fio) = "/^o^o-^co;/?) and the highest weight vector |Ao >= \P,P,P,0 > satisfies 

r/o|Ao >= 0. 



So we make the following ansatz ||T6 



\0;r>B 



Fn(r) 



|Ao >, 



E 



m 



2 — , imi 

32 7 2 



-a. 



m=l 



(3.20) 



(3.21) 



where /^m' /^m' /^m functions of the boundary parameter r. 

We can check that c^"*^^) plays a role of the Bogoliubov transformation 
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h*^ 4- ry h*^ 4- 



\m\ 



m 



-Fo(r);^*2 Fo(r) ^ /i*2 , ^*2 , M!^2 

m 



p-Fo{r)„ pFo{r) 



r'mi 



Keeping ( p.l8| ) in mind and following [0, |T9|, we find that the coefficients am, /^m? 



are 



/5^ = o, 



= ri^'" + ^^ 



-,2m 



\m\ 



(3.22) 
(3.23) 

(3.24) 



where the function 9m is defined by 

Om. = \ 



1 if m is even 
if m is odd 



Moreover following [|T^] one can check that ?7o|0 >b= 0, namely, the boundary state |0 >b 
G V{fio), as required. In the derivation, the following relation are useful 



ght+{rV;-|)|0;^>^ = e'^r(?5^;-i)|0;r >B, 



|0 >B = (l-cur)-'e 



10;r >B, 



|0;r>B. 
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Similarly, the dual state ^ < r; 0| G V*{fio) can be constructed 



B<r;0\ 
Goir) 



1 oo 

r) ^ y r ]2 \"'m"'m ^ "'m"'ni ' ^m^m ( 
oo 



m=l 



where 



^2 



_m — — 



q 2 + q- 



\m\ 



\m\ 



3.2 The general boundary states 

Noting that the boundary K-matrix K{z) have the following properties 

/ 1 







V 







(3.25) 



(3.26) 



(3.27) 



we may define | — l;r >b= $i(r~-'^)|0; r >b \r .rg-^- One can check that such | — l;r >b 

satisfies (|3.18|) with a = —1. Recursively, we can construct the general boundary state 
\a]r >B from \0]r >b by the following recursive relations. 



rq^ >B= $i(r ^)|a + 1; r >b, 



\a: r 



>B= q^^''<^Ur-'q^)\a - l;rq^ >b • (3.28) 



We have used the second invertibility relation ( [A.16|) . Similarly, we can obtain the dual 
boundary states ^ < r; a| from ^ < r; 0| by the recursive relations. 



B < r; a|<l>^(r) = b < rq'^; a — 1|, b < a\ = b < rq'^; a — l|$i(rg ^)g" 



-2pi 



(3.29) 



4 Correlation functions 

The aim of this section is to calculate the one-point functions < E^^j >q,: 

(1) _ B < r;a\EQ\a;r >b 



" B < r; a|a; r >B 

The generalization to the calculation of multi-point functions is straightforward. Thanks 
to the recursive relations (|3.28| ) and (p.29|) , it is sufficient to calculate < E^^j >o. Thus 

in the following we restrict ourselves to the calculation of < E^^j >q. 
Define 

<l dz f{z) = /_!, for formal series function f{z) = ^ fnz"'- 
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By the bosonic realization of Drinfeld currents of [/q[sZ(2|l)], ( |A.7|) -( |A.10[ ) and the normal 
ordering relations in the appendix A, we obtain the integral expression of the vertex 
operators M 



0lW 



: e 



-c(tog;0) 



+ 



jg-^2('?^^;~^)"'^2(l«;-|) + c(lj22;0)gi7rhj 



^(l_Hi£)(l_H2) 
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^hl{qz:-\) . gJTrag 



dw- 
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/ dwi i dw- ^ 
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w 

-c{wiq-fi) _ -c(«)ig-l;0) i , 

X : Qhl(qz--^)-hx(w-^)~h2{wv~^)^i-Kaf^ . 



Since rio\0,r >= 0, one may set 

B < r;0|$*(zi)$j(z2)|0;r >b _ b < r;0\(f)*{zi)(f)j{z2)\0;r >b 



< r; 0|0; r >b 



< r; 0|0; r >b 



(4.1) 



then < >o= -(-l)b]p„(l,l). 

The bosonization formulae ( |A.7| )- (|A.1CI|) of the vertex operators immediately imply 

P.,(z„ = Z2) , B<r;O|0r(z,)0.(..)|O;r>. _ 

1, 2j tj 1, 2; tj ^ ^ r;0|0;r >B 

Using the technique in [0, |21|] (see equation ( |C.4| )), after tedious calculation, we get 
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where 



n\ " ml ' 



„2n^n „—2n^—n ^^n, ,n „n, ,—n 

^ ~ ^ [n] """^ [n] ^ [n] ^"""^ [n] ' 



oo oo 
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We now derive the difference equations satisfied by the one-point functions. By ( |2.11|) 
and ( ]A.15| )-( [A.16D , one obtains 



j 



(4.6) 
(4.7) 



By ( [A.14| ), one derives the exchange relations 



kl 



Z2 



where R (z) = R-^^'*^-'^{z). 

Using (gT|)-(gJ|), (PI)- 
tions 



(|A.14|) and ( |A.7] )-( |/O!0| ), we get the difference equa- 



F,{z,q-\z2) = E (-i)^'="'i+[^^+f^'^+t"ii^K^i^"')^(^r'^2-'?')S 

j,k,l,m,n 



xKr{z,q-')R{^)llF^iz,,Z2), (4.9) 
Z2 
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Fi{zi,Z2q'^) 



(4.10) 
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A Appendix A 

A.l Bosonization of [/g[s/(2|l)] 

In this section, we briefly review the bosonization of ?7q[s/(2|l)] at level-one and the 
corresponding vertex operators j^ , |^. 
The Cartan matrix of ?7g[s/(2|l)] is 



where i,j = 0, 1, 2. 

In terms of the Drinfeld generators: {d, X^'^, h\^^ (i^'*)^-^, 7^^/^|z = 1,2, m G Z,n G 
Z^o}) the defining relations of ?7g[s/(2|l)] read 



[X^l^,X^'^]^±a^j + [X„|^;^,X^'*]^±a,^ = 0, for aij ^ 0, 
where [m] = ^^f^, [X^Y]/: = XY -{-1)™^^^YX and [X,Y]i = [X,Y]; the Zs-grading 



of Drinfeld generators are : [^m'^] = 1 for m G Z and zero otherwise. 

Introduce the bosonic q-oscillators |^ {a^, a^, 6„, c„, Qa^,Qa2,Qb,Qc \n G Z}, which 
satisfy the commutation relations 




7 is central, [K\hi] = 0, [d,K'] = 0, [d,hl] 



mhi 





m 
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[?77.] ^ 

[bm,bn] = —Sm+nfi 5 [bo,Qb] = " 1, 

m 

Cm)C„ = Sm+n,0 ; Cq, = 1. 

m 

Define the generating functions for tlie Drinfeld basis by Xf(z) = Smez ^m'*-^"™"^? 
introduce by setting = q^o. Define Qhi = Qa^ — Qa'^, Qh^ = Qa? + Qb and hi{z] [3) 
by 

K{z- (3) = -Y. ^g-'^l"!^-" + + h}. In z, (A.l) 

wliere /3 is a parameter. Otlier bosonic fields are defined similarly. 

The Drinfeld generators at level-one are realized by the free boson fields as P3 



hi = alq-\-\'' - alq\-\l\ hi = a^g"!-!/^ + m G Z, 

X^{z) = ± : e^'^i^^'^^) : e±*™o, X+{z) =: e'^2{^;i)gc(^;0) . g-^aj^ 

where dzf{z) = ^^'^^^Zq-I^z ■ ^ ■ stands for the usual normal ordering of O. 

Consider the bosonic Fock spaces -Fai,A2,A3,A4 5 generated by a*_^, b^m, c^m ijn > 0) over 
the vacuum vectors |Ai, A2, A3, A4 >, 

-Fai,A2,A3,A4 = C[a*_i, a*_2, b-i, ..; c_i, ..]|Ai, A2, A3, A4 >, (A. 2) 

where 

a^|Ai, A2, A3, A4 >= 0, 6^|Ai, A2, A3, A4 >= 0, c„|Ai, A2, A3, A4 >= 0, for m > 0, 
|Ai, A2, A3, A4 >= e^^<3"i+^2^«2+^3Q6+A4Qc|o^ 0, 0, > . 

Introduce the following spaces 

= -F)3+i,/3-i+i,/3-a+j,-a+j- (A.3) 

It can be shown that the bosonized action of f/g[sZ(2|l)] on -F(q;/3) is closed. To obtain the 



irreducible subspaces in it convenient to introduce a pair of fermionic currents pq. 



rjiz) = E v^z—' =: e<'-''^ :, ^z) = ^ ^n^"" =: e" 



-c(z;0) . 



The mode expansion of t]{z),C,{z) is well defined on F(^a;i3) for a G Z, and it satisfies the 
following relation 

Since rjo commutes (or anticommutes) with [/q[s/(2|l)], rjo plays the role of screening 
charge and ?7o^o quahfy as the projector from F(^a;i3) to the Kernel of r/o- Set A^ = (1 — 
a)Ao + aA2, a G Z, where Ai{i = 0, 1, 2) are the fundamental weights of ?7g[s/(2|l)], and 

' A„, a = 0,1,2 
= < for a > 2 . (A. 4) 

Xa for a < 

Define V{^a) = Vo^oF(a,f3-a)- Following |^, V{fia) G Z) are the irreducible highest 
weight f/g[sZ(2|l)]-modules with the highest weight /Iq. 
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A. 2 Level-one Vertex operators 

Let V^(A) be the highest weight f/g[s/(2|l)]-module with the highest weight A. Consider 
the following intertwiners of f/g[s/(2|l)]-modules: 

They are intertwiners in the sense that for any x G f/g[s/(2|l)], 

e{z)-x = A{x)-e{z), e(z) = $(^),$*(2), (a.s) 

the grading of these operators is :[0(2;)] = 0. is called type I (dual) vertex operator 
0]. We expend the vertex operator as 

$(^) = ^ <l>{z)j^v, , <^*{z)= '^*{z)j^vf. (A.6) 

i=l,2,3 j=l,2,3 

Define the operators (j)j{z), (f)*{z),ipj{z) and ip*{z) (j = 1,2,3) by 

03(2) =: e-/i2(9^^;-i)W^;0) . g-i-ag^ ^^^7^ 

M^) = -[03(^),Xo-'\-i, M^) = [M^),XoX (A.8) 
0*(2) =: e'^i^''^'-^) : e^™', (A.9) 
= -g-M</>tW,^o"'']., 03(^) = q^'m{z),X^\, (A.IO) 

where /i;;^ = -/i^, /i^^ = -/i^ - ^/i^ and Q/,*! = -Q?,2, Qh*2 = -Q^i - 2Qh2. Since 
the operator (f)i{z), (f)*{z) commute (or anti-commute) with rjo, we define 

H^) = r]oiMz)mio, = r]oio(l)l{z)r]oio. (A.ll) 

According |2^, 0, the vertex operators ^{z) and ^*{z) ( |A.6D given by (|A.11|) are the only 
type I vertex operators of [/g[s/(2|l)] which intertwine the level-one irreducible highest 
weight f/g[s/(2|l)]-modules V{^ia) (« G Z) 

<l>(^) : V{^l^) K(/i,_i) ® V;, $*(^) : \/(/i„) — . l^(/i„+i) ® VT- 

It is shown [0] that the above vertex operators satisfy the graded Faddeev-Zamolodchikov 
algebra 



^,{z,)^lz,) =Y.R{^tMz,)^,{z,){-\t^^, (A.12) 

kl ^2 

<|.*(^2)$*(^i) = E^(-)fc^'^fe(^i)'^K^2)(-l)f^"^'l, (A.13) 



fci ^2 

$,(^2)$r(.i) = E^(-)g"^u^i)<^>/(^2)(-i)f'=iw, (A.14) 

feZ ^2 

where -R(-2) = -R~^''**^(-2). Moreover, the vertex operators having the following invertibility 
relations 

*.W*;iv(A.) = -(-l)[^']5.,^rfy(A„), (A.15) 

A: 

$*(^g2)<l>,(z)|v.(A„) = 5.,g^'^'^rf|y(A.), (A.16) 
Er''"=*.(^)*fc(^g')|y(A.) = ^rf|y(A„)- 
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B Appendix B 

In this appendix, we give the normal ordering relations of fundamental bosonic fields: 



hi{zi;l3i) hi{z2;l32) . 



Mzi-A) h2{z2-A) ^ I . Mzi-A) h2{z2;(32) . 

h2izi;Pi) hi{z2;P2) _ \ . „h2{zi;Pi) hi(z2;P2) . 

^h2(zi-A)^h2(z2-A) ^. gft2{2i;/3l)gh2(22;/32) . 

gfti(2i;/3i)gh*(22;/32) ^ _ q-(f^i+f^^) z^^iJ : g'^'^^i-'^i^e^*^^"'*^ :, 

g/i,*(2i;/3i)g/ij(z2;/32) _ _ ^-(/3i+/32)^^^5i3 . g/J* (2i;/3i)g/»j(22;/32) 

gft];(2i;/3i)gft5;(22;/32) ^. ^hl{zi;f3i) ^hi{z2-A) 

hl{zi;f3i) h*{z2;f32) _ ^ . ^hUziM h*{z2-A) . 

zi- g-(/3i+/32)^2 • ^ 

/i*(2i;ft) Ji*(22;/32) ^ 3^ . ph*{zv,Pi)phl{z2;P2) . 

z,- q-iP^+Mz^ • ^ 

phl{zi;l5i) hl{z2\l32) _ \ . hl{zi;0i) hl{z2;l32) . 

{zi-q-(^^+P^)+^Z2){z,-q-(^^+^^)-^Z2) 

gC{2i;/3l)gC(22;/32) _ _ ^-(/^l+/^2)^2^ . gC(2i;/3l)gC(22;/32) . _ 



C Appendix C 

We here summarize the formulas concerning coherent states of bosons which have been 

used in section 5. 

The coherent states \(^, C^, > and < C^, in the Fock space -F(o;/3) and its dual 
space -^(o-^) are defined by 

\C\ C^ >= exp ( ^ ^ T^Ch-m + E T^Clc_^| 1/5, /3, P, >, (C.l) 

lm=l i=l L^J m=l L"^J J 

lm=l i=l ['"■J m=l L'"! J 

where and = 1> 2, 3, m = 1, 2, • • ■) are complex conjugate parameters. 
Noting that 

hU(3,(3,(3,0>=0, <f3,f3,f3,0\h'_^ = 0, t = l,2, m > 1, 
Cm\(3,(3,(3,0>=0, < A/5,A0|c_^ = 0, m> 1, 

one can easily verify 

c^ic\ c^ c^' >= dic\ c^ >, < c\ c^ c^ic-„. = c;^ < c\ c^ c^i- 

One can also show that the coherent states {|C^,C^,C^ >} (resp. < C^,C^,C^|}) form a 
complete basis in Fock space -F(o;/3)) (resp. FZ^.p\)- Namely, one can verify the completeness 
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relation 



idi 



(0;/3) 



f TT ^Sm^^m 
J 11 [m]2 



[aijm][m] • 



m=i ^ 'iet(^ 

x|c^c^c^><c^c^c^ 



^ 2 x,.(,n)m 

m=l i,j=l 



\m\ 



(C.3) 



where Kij{n) is a 2 x 2 matrix satisfying 



Y^Kii{n)[aijn\ = 

1=1 



One may also derive the following identity 



n 



dCm^Cm ^Cm^Cm ^Cm^Cm. 
m V m / 



-j^ oo 

o ^ -^TO (Cn' Cm' Cm' Cm' Cm' Cm) "^^i 



m=l 



Sm 

V Cm y 



~ vSm' Sm' Sm' Sm' Sm' Sm/*-^' 

oo 



m=l 



mil [2^^ m 



i3 ^" i3 



m-^m '-^nJ f ' 



(C.4) 



where Am are invertible constant 6x6 matrices and Bm are constant 6 component vectors. 
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